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HIGHLIGHTS 


►  Nanostructured  thin-film  ternary  PtCoMn  cathode  catalyst. 

►  Pressure  and  relative  humidity  dependence  of  oxygen-reduction  reaction  kinetics. 

►  39.5  kj  mor1  activation  energy  for  ORR  kinetics. 

►  0.13-0.25  A  mgn  mass  activity  at  0.9  V,  1  atm  H2  and  02,  80  °C  and  100%  RH. 

►  1.7— 2.0  mA  cni],2  specific  activity  at  0.9  V,  1  atm  H2  and  02,  80  °C  and  100%  RH. 
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The  kinetics  of  the  oxygen  reduction  reaction  (ORR)  on  nanostructured  thin-film  (NSTF)  ternary  PtCoMn 
catalyst  was  investigated  in  50  cm2  proton-exchange  membrane  single  cell  fixtures.  The  investigation 
covered  a  wide  range  of  oxygen  partial  pressures  (0.2-0.5  atm),  temperatures  (30-90  °C),  relative 
humidities  (25—100%),  and  Pt  loadings  (0.05—0.2  mg  cm~2)  in  the  cathode  catalyst.  The  kinetic  data  were 
derived  from  cell  polarization  curves  measured  during  Galvanodynamic  scans  from  near  open-circuit  to 
a  high  current  density  and  back  to  near  open-circuit.  The  kinetic  data  were  correlated  with  a  single  Tafel 
equation  and  a  transfer  coefficient  that  is  a  function  of  the  relative  humidity  (RH).  The  correlation 
indicates  a  0.36  order  for  the  02  partial  pressure,  39.5  kj  mol  1  activation  energy  for  the  temperature 
dependence,  and  an  additional  0.9  order  for  RH  dependence.  The  calculated  mass  activities  and  specific 
activities  at  0.9  V  Ohmic  resistance  corrected  cell  voltage,  1  atm  H2  and  02  partial  pressures,  80  °C  and 
100%  RH  agree  well  with  the  measured  values  of  0.13-0.25  A  mgpt1  and  1.7-2.0  mA  cmP12  for  the  cells 
with  different  Pt  loading. 

Published  by  Elsevier  B.V. 


1.  Introduction 

Nanostructured  thin  film  electrocatalysts  have  shown  promise 
in  addressing  the  performance,  cost  and  durability  barriers  that 
must  be  overcome  for  mass  commercialization  of  automotive  fuel 
cells.  Several  studies  by  Debe  and  associates  [1-3]  have  demon¬ 
strated  that  the  NSTF  catalysts  can  have  significantly  higher  specific 
activity  and  durability  than  the  conventional  carbon  supported 
highly  dispersed  Pt  catalysts.  The  NSTF  catalysts  are  formed  by 
vacuum  sputter-deposition  of  catalyst  alloys  onto  a  supported 
monolayer  of  highly  oriented  crystalline  organic-pigment  whiskers 
[1],  The  support  whiskers  have  a  high  aspect  ratio  (20-50), 
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0.6—2  pm  length,  rectangular  lath-like  morphology,  55  ±  12  nm 
width  and  27  ±  7  nm  thickness,  and  3—5  x  109  cm-2  a  real  number 
density  [1  ].  The  whiskers  have  been  coated  with  pure  Pt  as  well  as 
binary  and  ternary  Pt  alloys  [2,3].  The  PtCoMn  ternary  alloy 
(nominal  68:29—30:2—3  atomic  ratio)  is  of  particular  interest  as  it 
has  been  evaluated  extensively  for  performance  and  durability 
under  potential  cycling  conditions  (cyclic  voltammetry  with  H 2/ 
N2),  long  holds  at  1.5  V  and  higher  potentials  (support  stability), 
and  startups  and  shutdowns  in  H2/air  [4,5], 

The  purpose  of  this  work  is  to  characterize  the  activity  of  the 
ternary  Pt0.6sCo0.3Mn0.02  NSTF  catalyst  for  the  oxygen  reduction 
reaction  (ORR)  in  H2/air  fuel  cells.  The  intent  is  to  determine  the 
ORR  activity  as  a  function  of  the  oxygen  partial  pressure,  temper¬ 
ature,  relative  humidity  and  Pt  loading.  Previous  studies  on  the  ORR 
activities  of  Pt  and  PtCoMn  catalysts  on  NSTF  whiskers  have  used 
the  rotating  ring-disk  electrode  (RRDE)  method  in  0.1  M  FICIO4  at 
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room  temperature  to  60  °C  [6],  These  studies  have  shown  that  the 
area  specific  activity  of  Pt/NSTF  is  similar  to  that  of  bulk  Pt  and  that 
the  activity  of  the  ternary  PtCoMn  catalyst  is  about  twice  the 
activity  of  Pt/NSTF.  The  observed  enhancement  in  the  ORR  activity 
in  the  RRDE  measurements  is  similar  to  the  results  obtained  in 
50  cm2  H2/air  fuel  cells  [6-10], 

Our  approach  to  determine  the  ORR  kinetics  for  PtCoMn/NSTF 
from  the  polarization  curves  for  H2/air  fuel  cells  benefits  greatly 
from  similar  work  described  in  literature  on  dispersed  Pt/C  elec¬ 
trodes.  The  past  work  has  relied  on  holding  the  cell  at  a  constant 
potential  for  10  min  or  longer  to  ensure  that  the  catalyst  reaches  the 
corresponding  equilibrium  Pt  oxide  coverage  [11  ].  The  advantage  of 
using  long  hold  times  is  that  the  data  are  quite  reproducible;  the 
disadvantage  is  that  the  data  cannot  be  directly  used  to  analyze  the 
polarization  curves  that  are  generally  obtained  at  faster  scan  rates. 
We  used  a  slightly  different  approach  of  holding  for  2  min  at  each 
data  point  and  averaging  the  potentials  measured  during  down- 
scans  from  0.02  to  2  A  cm-2  and  reverse  upscans  from  2  to 
0.02  A  cm-2  (10  points  per  decade)  to  estimate  cell  voltages  at 
conditions  approaching  the  equilibrium  oxide  coverages.  Future 
work  will  explicitly  determine  the  oxide  coverage  in  the  ternary 
catalyst  as  a  function  of  the  cathode  potentials,  the  effect  of  the 
oxide  coverage  on  the  ORR  activity,  and  the  dynamics  of  Pt  oxida¬ 
tion  and  Pt  oxide  reduction. 


2.  Experimental 

Eight  50  cm2  active  area  single  cells  were  assembled  with 
a  24  gm  membrane  (850  equivalent  weight),  the  ternary  NSTF 
catalyst,  and  3  M  gas  diffusion  layers  (GDL)  into  a  Fuel  Cell  Tech¬ 
nologies  test  cell  containing  quad  serpentine  flow  fields.  The  GDL 
consisted  of  a  backing  paper  to  which  was  applied  a  hydrophobic 
treatment  and  an  MPL  (micro-porous  layer).  All  cells  had  a  Pt 
loading  of  0.05  mg  cm-2  in  the  anode.  Two  of  the  eight  cells  had  a  Pt 
loading  of  0.103  mg  cm  2  in  the  cathode.  The  Pt  loading  in  the 
cathode  in  the  other  cells  (two  each)  was  0.054,  0.146  and 
0.186  mg  cm-2.  All  cells  were  first  conditioned  using  a  “thermal 
cycling”  process,  described  in  detail  in  Steinbach  et  al.  [12],  which 
consisted  of  repeated  temperature  and  voltage  cycles  over  a  period 
of  2—3  days  until  stable  performance  was  reached.  Next,  oxygen 
reduction  activity  of  the  cathode  catalyst,  short  resistance  and  H2 
crossover,  and  electrochemically  active  surface  area  (ECSA)  of 
cathode  catalyst  were  measured,  after  which  a  series  of  tests  were 
conducted,  as  described  below. 

Table  1  lists  the  baseline  Pt  loading  in  the  electrodes  and  the 
operating  conditions.  One  of  the  cells  with  0.103  mg  cm-2  Pt 
loading  in  the  cathode  was  considered  as  the  reference.  The  refer¬ 
ence  operating  conditions  were  80  °C  cell  temperature,  1.5  atm  H2 
and  air  inlet  pressures,  65  °C  inlet  dew  point  temperatures,  and  H2 
and  02  stoichiometries  (SR)  of  2. 

Table  2  summarizes  the  test  matrix  with  eight  series  of  tests. 
Series  1  tests  varied  the  cell  temperature  from  75  °C  to  90  °C 
(P  =  1.5  atm,  SRc  =  2,  SRa  =  2)  with  the  inlet  dew  points  pre¬ 
determined  as  a  function  of  pressure,  temperature  and  anode/ 
cathode  stoichiometry  to  maintain  100%  relative  humidity  (RH)  at 
cell  exit.  In  all  tests,  humidification  water  was  injected  directly  into 


Reference  Pt  loading  in  electrodes  and  cell  operating  conditions.  The  outlet  RH  is 
100%  at  all  reference  operating  conditions. 
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the  anode  and  cathode  feed  streams  to  reach  the  set  dew  points. 
Also,  the  outlet  RH  was  not  measured  but  was  estimated  from  the 
cell  operating  conditions. 

The  backpressure  regulator  was  used  to  vary  the  cell  inlet 
pressure  in  Series  2  tests  from  1  to  2.5  atm  (T  =  80  °C,  SRc  =  2, 
SRa  —  2).  As  explained  later,  the  target  inlet  pressure  could  not  be 
reached  under  some  conditions  even  with  the  backpressure  valve 
completely  open  to  the  ambient.  As  in  Series  1  tests,  the  inlet  dew 
points  were  also  adjusted  to  maintain  100%  RH  at  cell  exit. 

Series  3  tests  varied  the  inlet  dew  points  from  dry  (no  water 
addition)  to  80  °C  (P  =  1.5  atm,  T=  80  °C,  SRc  =  2,  SRa  =  2).  Series  5 
tests  varied  the  02  stoichiometery  from  1.5  to  5  with  the  inlet  dew 
points  adjusted  to  maintain  100%  RH  at  cell  exit  (P  =  1.5  atm, 
T=  80  °C,  SRa  =  2).  Series  6  tests  varied  the  H2  stoichiometry  from 
1.2  to  5  with  the  dew  points  held  at  the  reference  value  of  65  °C 
(P  =  1.5  atm,  T  —  80  °C,  SRc  =  2).  Series  7  tests  were  designed  to 
investigate  the  cell  operation  during  warm-up  at  30,  45  and  60  °C 
with  dry  feeds  (P  =  1.5  atm,  SRc  =  2,  SRa  =  2).  Similarly,  Series  8 
tests  were  designed  to  investigate  cell  operation  during  idling  at 
low  pressures  and  high  02  stoichiometries. 

Series  4  tests  were  conducted  on  cells  with  different  Pt  loadings 
in  the  cathode.  These  cells  were  only  operated  at  the  reference 
conditions. 

To  check  for  cell-to-cell  variability,  all  tests  were  repeated  on  the 
companion  cells  with  identical  Pt  loadings  in  the  cathode. 

Fig.  la— h  present  the  polarization  curves  for  the  eight  series  of 
tests.  These  were  obtained  during  Galvanodynamic  scans  (CDS)  at 
cell  current  densities  varying  from  0.02  to  2  A  cm  2.  Only  the  data 
up  to  1  A  cm-2  are  plotted  in  Fig.  1  since  significant  mass  transfer 
effects  are  present  at  higher  current  densities.  Our  procedure  was 
to  start  the  measurements  near  the  open  circuit  at  0.02  A  cm-2, 
then  increase  the  current  density  in  steps  until  either  2  A  cm-2  was 
reached  or  the  cell  voltage  decreased  to  0.40  V,  followed  by  an 
upscan  back  to  0.02  A  cm-2.  The  cell  voltages  were  measured  at  ten 
current  steps  per  decade  but  the  maximum  step  size  was  limited  to 
0.1  A  cnT2.  The  cell  was  held  for  120  s  at  each  current  step  and  the 
cell  voltage  was  recorded  every  5  s.  Also,  the  cell  resistance  was 
determined  every  5  s  by  on-line  AC  impedance  measurements 
between  10  and  1  kHz,  from  which  the  high  frequency  resistance 
(HFR)  was  estimated  by  fitting  to  a  Randle’s  cell  model.  The  cell 
voltages  presented  in  Fig.  1  are  the  measured  values  averaged  over 
the  last  30  s  at  the  set  current  densities. 

Prior  to  the  start  of  the  experiments,  the  electrochemical  surface 
areas  were  determined  by  cyclic  voltammetry  (CV)  between  0.085 
and  0.65  V  (vs.  RHE)  at  100  mV  s^1  scan  rate,  with  H2/N2  at  70  °C, 
70  °C  dew  points  and  ambient  pressure.  ECSA  values  were  deter¬ 
mined  by  averaging  the  oxidative  and  reductive  Hupd  peak  areas 
from  100  CVs  and  assuming  a  210  pC  cmpt2  charge-to-area 
conversion  factor.  The  measured  ECSA  values,  listed  in  Table  3, 
decrease  from  12.4  m2  g-1  at  0.054  mg  cm-2  Pt  loading  to 
~  7  m2  g-1  at  0.186  mg  cm-2  Pt  loading,  and  are  five  to  ten  times 
smaller  than  the  ECSA  of  highly  dispersed  Pt/C  (55-70  m2  g-1  for 
Pt/C  with  0.4  mg  cm-2  Pt  loading). 

Table  3  includes  the  surface  enhancement  factors  (SEF)  obtained 
by  multiplying  the  measured  ECSA  by  the  Pt  loading.  The  SEF 
(cmpt.cm-2)  grows  with  increase  in  Pt  loadings  at  a  rate  that  is  less 
than  linear.  Also  included  in  Table  3  are  the  thicknesses  of  the 
catalyst  layers  as  measured  from  SEM  cross-sectional  images  of  the 
NSTF  catalyst-coated  membranes  [10].  These  NSTF  catalyst  layers 
are  only  about  one-tenth  as  thick  as  the  conventional  Pt/C  catalysts 
with  similar  Pt  loadings. 

Hydrogen  crossover  current  density  and  cell  short  resistance 
were  determined  by  measuring  the  plateau  currents  between  0.4 
and  0.6  V  cell  potentials  with  saturated  H2/N2  at  1.5  atm  total 
pressure  and  80  °C  cell  temperature.  Table  2  shows  that  the  H2 
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Test  Series  1:  Effect  of  temperature 

Pt  loading:  0.1(c),  0.05(a)  mg.cm"2;  Pressure:1.5  : 

Run  1.1 

Temperature, °C  75 

Dew  point  T(c),°C  55 

Dew  point  T(a),°C  55 

Test  Series  2:  Effect  of  pressure 

Pt  loading:  0.1(c),  0.05(a)  mg.cm"2;  Temperature: 

Run  2.1 

Dew  point  T(c),°C  73 

Dew  point  T(a),°C  73 

Pressure,  atm  1 

Test  Series  3:  Effect  of  relative  humidity 

Pt  loading:  0.1(c),  0.05(a)  mg.cm'2;  Temperature: 

Dew  point  T(c),°C  0 

Dew  point  T(a),°C  0 

Test  Series  4:  Effect  of  platinum  loading 
Temperature:  80  °C;  Dew  point:  65  °C(c),  65  °C(a);  P 
Run  4.1 

Pt  loading(c),  mg.cm"2  0.05 

Pt  loading(a),  mg.cm"2  0.05 

Test  Series  5:  Effect  of  cathode  stoichiometry 
Pt  loading:  0.1(c),  0.05(a)  mg.cm"2;  Temperature:  8 
Run  5.1 

Dew  point  T(c),°C  60 

Dew  point  T(a),°C  60 

SR(c)  1.5 

SR(a)  2 

Test  Series  6:  Effect  of  anode  stoichiometry 

Pt  loading:  0.1(c),  0.05(a)  mg.cm"2;  Temperature:  8 

SR(c)  2 

SR(a)  1.2 

Test  Series  7:  Cold  temperatures 
Pt  loading:  0.1(c),  0.05(a)  mg.cm"2;  Pressure:1.5  atr 
Run  7.1 

Temperature, °C  30 

Dew  point  T(c),°C  0 

Dew  point  T(a),°C  0 

Test  Series  8:  Idle  operation 

Pt  loading:  0.1(c),  0.05(a)  mg.cm"2;  Temperature:  8 

Pressure,  atm  1.05 

SR(c)  10 

SR(a)  2 


crossover  current  density  of  the  membrane  at  these  pressure  and 
temperature  conditions  was  about  2.4— 3.5  mA  cm-2  and  the  short 
resistance  of  the  eight  cells  varied  between  182  and  933  Q  cm2. 

Table  3  lists  the  mass  activities  of  Pt  in  each  cell  measured  in  H2/ 
02  at  80  °C,  1  atm  reactant  H2  and  02  pressures  and  100%  RH.  The 
mass  activity  measurement  was  made  by  first  reducing  the  cathode 
with  a  15  min  hold  under  H2/N2  at  ~  0.10  V,  after  which  the  cathode 
was  switched  to  02  and  the  cathode  potential  was  then  immedi¬ 
ately  held  at  0.90  V  for  20  min.  The  listed  mass  activity  values  are 
the  current  densities,  adjusted  for  H2  crossover  and  shorting 
currents,  at  a  cell  voltage  of  0.9  V  after  17.5  min.  The  measured  mass 
activity  is  0.18  A.  mgpt1  for  the  reference  cell  with  0.103  mg  cm-2  Pt 
loading  and  is  somewhat  higher  in  cells  with  lower  Pt  loading.  This 
value  may  be  compared  with  the  0.13  A  mgpt1  mass  activity  of  47% 
Pt/C  (TKK)  reported  by  Neyerlin  et  al.  [11]  and  0.1  A  mgpt1  mass 
activity  of  TKK  PtCo  dispersed  on  corrosion  resistant  (CR)  carbon  as 
reported  by  Wagner  et  al.  [13],  Thus,  as  noted  in  earlier  publica¬ 
tions,  the  mass  activity  of  ternary  PtCoMn/NSTF  is  40%  higher  than 
Pt/C  in  spite  of  the  much  lower  ECSA  and  80%  higher  than  PtCo/CR- 
carbon  [6,7],  The  reason  for  this  result  is  the  much  higher  specific 
activity  of  PtCoMn/NSTF.  Table  3  shows  that  the  ternary  NSTF 
catalyst  has  a  specific  activity  between  1.7  and  2  mA.cmp2,  which  is 


6—8  times  the  reported  values  of  0.25  mA  cmp2  for  Pt/C  and 
0.3  mA  cmp2  for  PtCo/CR-carbon  [13]. 

3.  Kinetics  of  oxygen  reduction  reaction 

We  investigated  the  ORR  kinetics  by  analyzing  the  polarization 
data  for  low  current  densities,  i  <  0.5  A  cm-2,  where  the  mass 
transfer  effects  are  negligible.  We  averaged  the  cell  potentials 
measured  during  the  upscans  and  downscans  in  order  to  remove 
any  non-equilibrium  effects  caused  by  the  finite  Pt  oxidation 
kinetics.  We  assessed  the  effect  of  the  finite  kinetics  of  Pt  oxidation 
and  PtO  reduction  near  the  open-circuit  potentials  by  running 
polarization  tests  at  the  reference  conditions  with  10  min  hold  at 
each  data  point.  The  data  showed  only  small  and  random  differ¬ 
ences  between  the  cell  voltages  measured  with  2  min  and  10  min 
hold  times.  We  excluded  i  <  0.1  A  cm-2  polarization  data  while 
determining  the  ORR  rate  constants  for  three  reasons.  First,  our 
analysis  procedure  involved  normalizing  the  error  (difference  in 
the  modeled  and  measured  current  densities)  by  the  measured 
current  density,  which  skewed  the  results  by  placing  a  large  weight 
on  the  small  current  density  data.  Second,  the  flow  rates  of  feed 
gases  are  held  at  threshold  values  for  i  <  0.1  A  cm-2,  so  that  that  the 
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anode/cathode  stoichiometries  and  exit  dew  points  are  different 
from  those  for  i  >  0.1  A  cm-2.  Finally,  the  crossover  current  (ix)  can 
be  a  significant  fraction  of  the  total  current  for  i  <  0.05  A  cm-2  but 
we  did  not  specifically  measure  the  RH  dependence  of  ix.  We 
emphasize,  however,  that  whereas  the  i  <  0.1  A  cm-2  data  were  not 
included  in  determining  the  ORR  rate  constants,  they  were 
considered  in  the  subsequent  evaluation  of  the  kinetic  model. 

We  represented  the  cell  voltage  (£)  at  the  measured  current 
density  (i)  in  terms  of  the  Nernst  potential  (En)  calculated  at  the  cell 
temperature  (T)  and  the  partial  pressures  (P)  of  H2  and  O2,  the 
cathode  overpotential  for  ORR  (?j),  and  the  Ohmic  drop  due  to  HFR 
(Rq)- 

e  =  en  -  mQ  -  t] 


a  stepwise  procedure  for  determining  7,  io,  a,  and  /3.  We  first 
determined  7  by  using  the  average  polarization  data  in  Series  2  at 
80  °C,  50%  H2  and  02  utilizations,  and  100%  RH  at  cell  exit.  For  these 
conditions,  Fig.  2a  presents  the  ORR  overpotential  as  a  function  of 
the  current  density  for  different  air  inlet  pressures  in  the  Series  2 
tests.  Fig.  2b  indicates  that  the  data  in  Fig.  2a  are  well  correlated  if 
the  current  density  is  scaled  by  PqT  with  a  7  of  0.36.  Here,PX 
represents  the  average  partial  pressure  of  O2  in  the  channel,  and  is 
calculated  by  accounting  for  the  pressure  drop  ( AP)  in  the  flow  field 
and  the  O2  consumption.  The  following  correlation,  developed  from 
the  measured  inlet  pressure  in  tests  in  which  the  exit  pressure  was 
1  atm,  was  used  to  estimate  AP  (kPa)  in  terms  of  the  gas  volumetric 


[1] 

flow  rate  {V  in  m3  s  superscripts  a  and 
cathode,  respectively). 

c  denote  anode  and 

[2] 

A Pa  =  1.6  +  1.3  x  105V  +  1.2  x  109\72 

APC  =  1.8 +  2.3  x  105V  +  5.9  x  1091/2 

[5] 

where  Eq  is  the  standard  potential  for  the  fuel  cell  reaction 
assuming  that  the  product  water  is  liquid,  $  is  the  relative  humidity 
(1  for  100%  RH),  and  n  equals  2. 

E0  rn  1.23  +  0.000833(T  -  298)  [3] 

We  neglect  the  small  anode  overpotential  for  the  facile 
hydrogen  oxidation  reaction  [10]  and  assume  that  the  sluggish 
oxygen  reduction  reaction  on  cathode  follows  the  Tafel  equation 
with  a  single  Tafel  slope  over  the  cathode  potential  range  of 
interest. 

i  +  ix  =  i0APtLPtPl2  #exp  [4] 

where  ix  is  the  hydrogen  crossover  current  density,  io  is  the 
exchange  current  density  (A.cmpt2),  Apt  is  the  electrochemically 
active  surface  area  (ECSA,  cmpt.mg-1),  Lpt  is  the  Pt  loading 
(mg  cm-2),  7  is  the  order  with  respect  to  O2  partial  pressure,  (J>  is 
the  relative  humidity,  f)  is  the  order  with  respect  to  <t>,  and  a  is  the 
transfer  coefficient.  It  is  assumed  that  H2  that  crosses  over 
undergoes  H2  =  2H+  +  2e-  half— cell  reaction  on  the  cathode 
catalyst,  produces  a  mixed  cathode  potential,  and  leads  to  an 
internal  cathodic  current  density  (ix). 

3.1.  Cell-averaged  kinetics 

We  investigated  the  cell  ORR  kinetics  by  using  the  average  H2 
and  02  partial  pressures  across  the  cell  (see  Appendix  1 ).  We  found 
that  the  ORR  overpotentials  correlated  better  with  the  RH  at  cell 
exit  rather  than  the  inlet  or  average  RH  values.  We  formulated 


In  the  second  step,  we  investigated  the  temperature  depen¬ 
dence  of  io  by  analyzing  the  Series  1  data  at  75,  80,  85  and  90  °C, 
with  50%  H2  and  O2  utilizations  and  100%  RH  at  cell  exit,  see  Fig.  3a. 
We  included  the  Series  2  data  for  1, 1.25,  2  and  2.5  atm  inlet  pres¬ 
sures  by  using  the  7  values  determined  in  the  first  step.  We 
correlated  the  data,  plotted  as  semi-logs  in  Fig.  3b,  with  straight 
lines  and  estimated  a  from  the  slope  and  io  from  the  intercept  as  per 
the  following  equation  ($  =  1). 

V  =  §Fln[ii  +  k)P*}  -Uln(i°Aptipt<pP)  [6] 

Fig.  3b  presents  the  estimated  io  as  a  function  of  1  /T  and  leads  to 
the  following  Arrhenius  correlation  with  an  activation  energy  (E0.2) 
of  57.7  kj  mol-1  for  the  ORR. 

where  ior  is  the  exchange  current  density  at  the  reference 
temperature  (Tr  =  353  K). 

The  third  step  consisted  of  investigating  the  RH  dependence  of 
the  ORR.  Fig.  4a  presents  the  ORR  overpotentials  in  the  Series  3  tests 
for  different  exit  RH  at  80  °C  and  in  the  Series  2  tests  for  1, 1.25,  2 
and  2.5  atm  inlet  pressures.  As  per  Eq.  (6),  we  fitted  the  data  with 
straight  lines  on  semi-log  axes  and  used  the  slopes  to  determine 
a  as  a  function  of  0,  see  Fig.  4b. 

0.4665 

“  -  1  +  0.0926$  lS| 


Current  Density,  A-cnr* 


Scaled  Current  Density  ( i-P02  »“),  A  cm-* 


Fig.  2.  Derived  cell-averaged  ORR  overpotentials  in  series  2  data  (a)  and  determination  of  the  order  with  respect  to  P0l  (b). 
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Fig.  3.  Derived  cell-averaged  ORR  overpotentials  in  series  1  data  (a)  and  determination  of  the  activation  energy  E0l  for  temperature  dependence  (b) 


We  used  the  intercepts  in  Fig.  4a  to  determine  and  these  are 
presented  in  Fig.  4c.  The  power  law  fit  in  Fig.  4c  leads  to  the 
following  estimates  for  ior  (A.cmpt2)  and  /J. 


lor  =  5.7  x  10  ' 
0  =  1-3. 


[9] 


3.2.  Local  ORR  kinetics 


Whereas  the  above  correlation  for  the  average  ORR  kinetic  rate 
is  useful  for  characterizing  the  cell  performance,  it  is  desirable  to 
have  the  ability  to  estimate  the  ORR  overpotentials  for  the  local 
values  of  i,  T,  P0l  ,  and  <I>.  For  this  purpose,  we  developed  a  simple 
one-dimensional  model  to  estimate  the  local  values  of  i,P0l  ,  PHl 
and  <f>  from  the  measured  cell  polarization  curves.  The  model  solves 
the  following  set  of  equations  for  anode  and  cathode  pressures  and 
for  molar  flow  rates  (N)  of  H2,  02  and  H20  (subscript  w)  in  the 
anode  (superscript  a)  and  cathode  (superscript  c)  channels 
considering  the  source  terms  for  the  electrochemical  reactions  and 
water  transport  flux  (/™)  across  the  membrane  (superscript  m)  by 
electro-osmotic  drag  and  diffusion  (normal  direction  y). 


dPc  _  APC 
dx  ~  if 


dpa  _  Apa 

dx  I 


[10] 

[11] 

[12] 

[13] 

[14] 

[15] 


Fig.  4.  Derived  cell-averaged  ORR  overpotentials  ir 


i  (a)  and  determination  of  the  transfer  coefficient,  a  (b)  and  the  order  with  respect  to  RH,  fl  (c). 
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The  following  auxiliary  equations  define  and  relate?^,  Ph2 
and  0  to  the  molar  flow  rates.  Also,  the  current  density  i  is  given  by 
Eq.  (4). 


[16] 

[17] 

[18] 

_c  AL,  PL 
$  =  =  -jf 

Nws  Ps 

[19] 

Ala  Ps 

[20] 

PH  = 

2  £Ni 

[21] 

Po  = 

2  £Ni 

[22] 

The  following  equation  is  used  to  calculate 
a  function  of  0  (<I>C  of  4a). 

A  (Ac  or 

Aa)  as 

f  23.8204  -  28.5303  +  8.3602  +  5.92<Z>+ 0.1 
■*=  <  9.67  +  2.8(0-1) 

( 15.27 

(0<1) 

(0<3) 

(0>3) 

[23] 

In  the  above  equations,  H2  flows  in  the  anode  channel  along  the 
x  direction  (L  is  the  total  length)  whereas  air  flows  in  the  cathode 
channel  along  the  -x  direction,  superscript  s  represents  saturation, 
A  is  the  water  uptake  in  the  membrane,  pm  (2000  kg  m  3)  and  EW 
(850)  are  the  membrane  density  and  equivalent  weight,  and  Wm  is 
the  membrane  width.  Note  that  Ps  is  the  saturation  pressure  of 
water  vapor  and  is  a  function  of  temperature  only.  Similarly,  Nws 
and  Nws  are  the  molar  flow  rates  of  saturated  water  vapor  and  are 
related  to  Ps  through  Eqs.  (17)  and  (18). 

For  expediency,  we  assumed  that  the  electro-osmotic  drag 
coefficient  (pj  =  1)  and  the  water  diffusion  coefficient  (D 
l  =  1.5  x  10-6  cm2  s-1)  are  constant,  so  that  Eq.  (16)  can  be  inte¬ 
grated  analytically  to  obtain  the  A  profile  and  hence  the  water 
transport  flux  J™  across  the  membrane. 

1241 

where  5m  is  the  membrane  thickness. 

A  simple  first-order  finite-difference  scheme  was  used  to  inte¬ 
grate  Eqs.  (10)— (15)  simultaneously.  We  marched  from  anode  inlet 
(x  =  0)  to  anode  outlet,  i.e.,  the  equations  were  integrated  in  the 
forward  flow  direction  for  the  anode  stream  and  in  the  backward 
flow  direction  for  the  cathode  stream.  The  required  boundary 
conditions  for  the  anode  stream  were  determined  from  the  known 
molar  flow  rates  of  H2  and  H2O  in  the  anode  feed  gas.  The  required 
boundary  conditions  for  the  cathode  stream  are  unknown  and  had 
to  be  determined  iteratively.  Our  approach  was  to  use  a  hybrid 
method,  combining  steepest  descent  with  quasi-Newton  update 


technique,  to  determinelV02  andNw  at  x  =  0  such  thatN02  andNw 
calculated  by  solving  Eqs.  (10)— (15)  at  cathode  inlet  match  the 
experimental  values  (control  variables). 

We  used  an  optimizer  to  determine  the  combination  ofior,£o2. 
and  p  that  minimizes  the  root  mean  square  difference  between  the 
measured  and  modeled  current  densities  normalized  by  the 
measured  current  density  in  Series  1-3  data.  The  parameters  were 
allowed  to  vary  over  a  wide  range:  5  x  10_8-5  x  10-6  A  cmpt2  for  ior, 
20-80  kj  mol'1  for£02,  and  0.25-2.5  for  p.  Attempts  to  vary 
a  proved  unsuccessful  because  of  convergence  problems  arising 
from  the  extreme  sensitivity  of  the  calculated  overpotentials  to  a.  It 
was  also  not  possible  to  vary  y  because,  over  the  range  of  pressures 
for  which  the  data  were  available,  the  calculated  overpotentials 
were  far  less  sensitive  to  y  than  to  other  parameters.  For  these 
reasons,  we  retained  Eq.  (8)  for  a  and  kept  7  at  0.36.  We  used 
a  sequential  quadratic  programming  technique  for  solving 
nonlinear  constrained  optimization  problems.  The  following  were 
the  initial  values  of  the  various  parameters  appearing  in  Eq.  (4)  as 
determined  by  the  optimizer:  ior  -  6.8  x  10-7  A  cmpt2,  E02 
=  39.5  kj  mol-1,  and  p  —  0.9.  In  a  follow-on  attempt,  P  was  allowed 
to  vary  linearly  with  0  since  subsequent  analysis  suggested 
a  stronger  dependence  of  ORR  kinetics  on  RH  under  dry  conditions 
than  can  be  explained  with  a  constant  p.  The  optimized  value  of  p  is 
given  by  the  following  linear  function. 


/  0.9  +  1.5(0.65-®)  0  <  0.65 
\  0.9  0  >  0.65 


[25] 


For  illustrative  purposes,  Appendix  2  shows  the  calculated 
profiles  of  i,Po2  ,  Ph2  and  0  in  Run  1.2  at  0.4  A  cm-2. 

Fig.  5  compares  the  measured  and  modeled  cell  current  density 
with  the  above  parameters.  It  includes  data  from  the  tests  in  Series 
1,  2  and  3  with  different  T,  P  and  RH.  All  these  tests  were  run  with 
50%  H2  and  O2  utilizations  and  a  Pt  loading  of  0.103  mg  cm-2  in  the 
cathode  and  0.05  mg  cm-2  in  the  anode.  For  the  data  included  in 
Fig.  5,  the  standard  deviation  (s)  between  the  modeled  (im)  and 
measured  (L)  current  densities  is  22.1  mA  cm  2  and  the  correlation 
coefficient  ( R 2)  is  0.993. 

Even  though  a  direct  comparison  between  the  ORR  kinetic 
parameters  derived  here  for  the  ternary  PtCoMn/NSTF  with  those 
reported  in  the  literature  for  Pt/C  is  not  meaningful  because  of  the 
different  catalysts  (alloy  vs.  monometallic,  polycrystalline  thin-film 


-  i0r =6.8x1 0-7  Acrr>t-2 
:  y=0.36 
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Fig.  5.  Test  of  the  local  ORR  kinetic  correlation  against  data  from  series  1,  2  and  3. 
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vs.  crystalline  nanoparticles,  nm  vs.  pm  curvature,  etc)  and  supports 
(organic  vs.  carbon,  whiskers  vs.  particles,  etc),  it  may  be  useful  to 
discuss  the  associated  trends.  The  PtCoMn/NSTF  system  shows 
a  small  RH  dependence  of  the  Tafel  slope  ( 2.303RTIanF)  that  at  80  °C 
decreases  from  82  mV/decade  for  100%  RH  to  78.5  mV/decade  for 
50%  RH.  The  Pt/C  system  has  a  constant  Tafel  slope,  65  ±  2  mV/ 
decade  [11,14],  although  dependence  of  a  on  RH  has  been  observed 
at  120  °C  [15].  The  pre-exponential  term  in  Eq.  (4),  however,  varies 
as  4>0'9  for  PtCoMn/NSTF,  whereas  Pt/C  shows  no  dependence  on  RH 
for  0  >  0.75  but  significant  reduction  for  $  <  50—60%  [16]  possibly 
due  to  its  effect  on  the  H+  concentration  in  the  ionomer  under  very 
dry  conditions  (A  <  2—4).  The  PtCoMn/NSTF  system  has  no  ionomer 
in  the  catalyst  layer  other  than  the  membrane  material  that 
partially  fills  the  voids  between  the  catalyst-coated  whiskers  when 
the  CCM  is  formed.  Also,  the  ORR  kinetics  on  PtCoMn/NSTF  shows 
smaller  dependence  on  Pq2  (y  of  0.36  vs.  0.54)  and  temperature 
(activation  energy  of  39.5  kj/mol  vs.  67  kj/mol)  than  Pt/C  [17,18], 
possibly  implying  different  rate-controlling  steps  in  the  overall 
mechanism  [10], 

The  graphs  in  Fig.  6  show  the  validity  of  the  kinetic  parameters 
with  respect  to  single  variables,  i.e.,  temperature  (and  RH)  in  Fig.  6a 
(Series  1  data),  pressure  (and  RH)  in  Fig.  6b  (Series  2  data),  and  RH 
in  Fig.  6c  (Series  3  data).  Fig.  6  shows  reasonable  match  between 
the  model  and  data  if  ior  is  slightly  adjusted  (by  -4%  to  6%)  to 
account  for  variations  in  measured  polarizations  in  the  different 
series  of  tests  (see  Fig.  11).  The  standard  deviation  and  R2  are 


13.1  mA  cm-2  and  0.9977  for  data  with  varying  temperatures, 
22.5  mA  cur2  and  0.9928  for  data  with  varying  pressures,  and 
22  mA  cm  2  and  0.996  for  data  with  varying  relative  humidities. 

Fig.  7  tests  the  applicability  of  the  model  parameters  in 
capturing  the  effects  of  cathode  and  anode  stoichiometry.  It 
includes  data  from  the  tests  in  Series  5  and  6,  which  were  not  used 
in  determining  the  optimum  choice  of  the  ORR  kinetic  parameters. 
We  found  that  ior  has  to  be  increased  by  ~  15%  at  SRC  -  5  and  ~40% 
at  SRc  =  10  (Fig.  8b)  in  order  to  get  a  good  match  between  the  model 
and  the  data  (s  =  34  mA  cm-2,  R2  =  0.9836).  This  increase  in  ior  at 
high  SRc  may  be  related  to  the  resultant  high  pressure  drop  that 
enhances  air  flow  under  the  lands  in  the  serpentine  flow  field.  It  has 
been  suggested  that  this  enhanced  air  flow  assists  in  removing  the 
liquid  water  from  the  GDL  under  the  land,  relieves  liquid  saturation 
in  the  cathode  catalyst  layer,  and  thereby  improves  the  catalyst 
activity  [19].  We  determined  an  empirical  factor  (/)  that  describes 
the  enhancement  in  ior  as  a  function  of  the  cell  cathode  stoichi¬ 
ometry  and  the  local  current  density. 


/  =  l+(fr(SR)-l)/2(i) 

MSR)  = 


/2(0={1 

lo. 


1.0  +  0.0033(SRC  —  2)z 
i>  0.125 


-0.0233(SRc-2) 


SR c>10 
SRC>1.75 


[26] 
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Fig.  6.  Test  of  the  local  ORR  kinetic  correlation  against 


data  with  varying  temperature  (a),  series  2  data  with  varying  pressure  (b)  and  series  3  data  with  varying  RH  (c). 
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Fig.  7.  Test  of  the  local  ORR  kinetic  correlation  against  series  5  data  with  varying  cathode  stoichiometry  (a),  and  series  6  data  with  varying  anode  stoichiometry  (b). 
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In  our  interpretation,  whereas  the  kinetic  parameters  a,  (?,  y,  ior 
and  EQl  are  unique  to  the  PtCoMn/NSTF  catalyst,  the  factor  /  is 
specific  to  the  flow  field  and  is  likely  to  change  with  a  different  cell 
design. 

As  expected,  the  match  between  the  model  and  data  is  better  in 
Fig.  7b  (s  =  9  mA  cm-2,  R2  —  0.9989)  than  in  Fig.  7a  since  the  anode 


stoichiometry  influences  the  ORR  kinetics  only  indirectly,  i.e., 
through  water  balance. 

Fig.  8a  provides  a  test  of  the  applicability  of  the  ORR  kinetic 
parameters  at  low  temperatures  that  the  cell  may  experience 
during  startup  (Series  7  data)  with  dry  feeds.  The  comparison  is 
restricted  to  0.1  A.cm-2  at  30  °C,  0.15  Axm-2  at  45  °C,  and  0.25 
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A.cm  2  at  60  °C  since  the  mass  transfer  effects  are  evident  at  higher 
current  densities  at  these  temperatures.  The  data  in  Fig.  8a  suggest 
that  an  activation  energy  of  39.5  kj  mol-1  adequately  reflects  the 
effect  of  temperature  on  ORR  kinetics  over  the  range  of  30—90  °C 
(s  =  2.3  mA  cm"2,  R2  =  0.9959). 

Fig.  8b  measures  the  “goodness”  of  the  correlation  for  low 
pressure  and  high-SR  data  (Series  8  for  idling  conditions,  100%  inlet 
RH).  We  see  that  the  empirically  derived  /  function  adequately 
describes  the  enhancement  in  ior  at  high  stoichiometries. 


3.3.  Discussion  of  the  results 

Fig.  9a  compares  the  mass  activities  measured  during  Series  4 
tests  with  the  current  densities  calculated  from  our  model  for  0.9  V 
IR-corrected  cell  voltage,  P0l  =  1  atm,  and  <i>  -  1.  As  stated  earlier, 
the  measured  mass  activities  correspond  to  current  densities 
measured  20  min  after  holding  the  cell  at  0.9  V.  Fig.  9b  gives 
a  similar  comparison  between  the  modeled  and  measured  specific 
current  densities.  Fig.  9a  and  b  include  data  from  cells  (pairs)  with 
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Cell  Voltage,  V 


Fig.  11.  Variability  in  data  from  six  repeats  of  tests  in  series  1,  2,  3,  4,  5  and  6  at  reference  conditions.  Ail  data  from  the  reference  cell  with  0.103  mg  cm  2  Pt  in  the  cathode. 


different  Pt  loadings.  The  favorable  comparisons  suggest  that  the 
ORR  kinetic  parameters  (a,  ff  and  y)  can  be  used  for  PtCoMn/NSTF 
cathode  catalysts  with  different  Pt  loadings  if  the  variation  in  ECSA 
with  Pt  loading  is  taken  into  account. 

Fig.  10  compares  the  modeled  and  measured  polarization  curves 
for  the  seven  series  of  tests  conducted  on  the  reference  cell.  There  is 
generally  a  good  agreement  between  the  modeled  and  measured 
cell  voltages  for  current  densities  between  0.1  and  0.4  A.citT2 
except  under  very  dry  conditions,  such  as  at  25%  RH  in  Series  3  data, 
where  there  is  additional  voltage  loss  due  to  membrane  dry-out. 
The  discrepancy  at  higher  current  densities  is  attributed  to  mass 
transfer  effects  that  will  be  analyzed  in  a  follow-on  paper.  The  mass 
transfer  effects  (possible  catalyst  layer  flooding)  are  evident  even  at 
0.1  A.cnrr2  for  cell  temperatures  below  45  °C  (Series  7  data).  The 
model  correctly  calculates  the  ordering  of  cell  voltages  with 
temperature  (£(75  °C)  >  £(80-85  °C)  >  £(90  °C))  in  Series  2  data. 
The  model  also  calculates  the  increase  in  cell  voltage  with  increase 
in  SRc  for  current  densities  between  0.1  and  0.4  A  cm-2,  as 
measured  in  Series  5  data,  and  the  observed  reordering  in  cell 
voltages  with  respect  to  SRc  at  lower  current  densities. 

It  is  worth  mentioning  that  the  inlet  pressure  is  the  control 
variable  except  at  low  pressures  (1—1.05  atm  data)  and  high  current 
densities  or  high  cathode  stoichiometries.  In  these  instances,  the 
outlet  pressure  is  1  atm  and  the  inlet  pressure  is  raised  to  produce 
the  specified  flow  rate,  see  Eq.  (5).  Thus,  the  cell  inlet  pressure  for 
the  Series  8  data  identified  as  1.05  atm  and  SRc  =  10  can  be  higher 
than  for  the  1.2  and  1.3  atm,  lower  SRc  data,  and  may  show  higher 
cell  voltages  at  current  densities  >0.1  A.cm-2. 

The  cell  voltages  calculated  from  our  correlation  are  within 
5—10  mV  of  the  measured  values  for  current  densities  lower  than 
0.1  A  cm-2.  In  taking  the  polarization  data  at  current  densities 
between  0.02  and  0.1  A  cm-2,  the  H2  and  air  flow  rates  were  held  at 
the  set  threshold  values,  100  slpm  for  H2  and  185  slpm  for  air.  Cell 
operation  at  high  air  stoichiometry  results  in  the  relative  humidity 
staying  close  to  the  inlet  value,  which  can  be  as  low  as  40%  in  the 
Series  2  test  at  75  °C  and  the  Series  5  test  at  1.5  nominal  SRc  (see 
Table  2).  We  found  much  better  agreement  between  the  model  and 
low-RH  data  for  <0.1  A.cm-2  current  densities  if  /?,  the  order  with 
respect  to  RH,  is  allowed  to  increase  for  4>  <  0.65,  as  in  Eq.  (25). 

The  polarization  data  for  the  reference  conditions  in  Table  1 
were  taken  at  different  times  in  Series  1,  2,  3,  5  and  6  tests.  We 
have  used  these  data  to  establish  the  variability  in  our  measure¬ 
ments.  Fig.  11a  indicates  that  the  measured  cell  voltages  at  the 
reference  conditions  are  reproducible  within  ±5  mV  and  have 
a  standard  deviation  of  2.1  mV  over  the  current  density  range  of 
0.02—0.5  A  cm-2.  Fig.  lib  presents  the  corresponding  variation  in 
the  current  density  estimated  by  using  a  Tafel  slope  of  82  mV/ 
decade  in  Eq.  (4).  The  standard  deviation  for  the  estimated  current 


densities  is  18.2  mA  cm-2  over  a  range  of  0.73-0.87  V  IR-corrected 
cell  voltages.  Our  kinetic  correlation  has  a  comparable  standard 
deviation  of  9—22.5  mA  cm-2  for  all  the  data  taken  in  Series  1, 2, 3, 5 
and  6  tests  along  with  ±6%  variation  in  ior.  Thus,  the  “goodness"  of 
the  ORR  correlation  is  consistent  with  the  accuracy  and  reproduc¬ 
ibility  of  the  polarization  data. 

Further  work  is  required  to  verify  and  improve  the  kinetic 
model.  The  model  reliability  will  improve  with  more  accurate  and 
reproducible  data  (see  Fig.  11).  Additional  data  are  needed  to  test 
the  RH  dependence  of  ORR  kinetics  (a  and  (?)  and  the  effect  of  high 
SRC  on  i0. 

4.  Summary  and  conclusions 

The  kinetics  of  the  oxygen  reduction  reaction  (ORR)  on  PtCoMn/ 
NSTF  can  be  correlated  with  the  following  simple  Tafel  equation 
with  a  single  value  of  the  transfer  coefficient  (a). 

i  +  ix  =  i0AnLPtPl2  4>pexp 

io  =  ior/ (SRc)  exp  [  -  %•  (f  -  ] . 

i0r  =  6.8x10  7A.cmpt2,  Eq2  =  39.5kJ.mol1,  Tr  =  353K 

y  =  0.36 

a  f  0.9  +  1.5(0.65  -  <2>)  <5  <  0.65 
P  /  0.9  <P  >  0.65 

0.4665 

“  -  1  +  0.0926$ 

/  =  1.  +  (ft (SR)  -  l)/2(i) 


/i(SR)  » 


1.4  SRC>  10 

1 .0  +  0.0033(SRc  -  2)2+0.0233(SRc  -  2)  SRC>1 .75 


/2(0 


i>0.125 
i<  0.125 


We  regard  the  kinetic  parameters  ior,  ol,  /?,  y  and  Eo2  as  specific  to 
the  PtCoMn/NSTF  catalyst  tested  and  the  factor  /  that  accounts  for 
the  increase  in  catalyst  activity  at  high  stoichiometries  as  an  attribute 
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of  the  flowfield  design.  The  ECSA  (Apt)  is  a  function  of  the  Pt  loading 
in  the  cathode  catalyst,  decreasing  from  12.5  ±  0.1  m2  g-1  for 
0.054  mg  cm-2  Pt  loading  to  7.11  ±  0.1  m2  g_1  for  0.186  mg  cm-2  Pt 
loading.  The  order  (/?)  for  RH  dependence  increases  sharply  for  $ 
below  0.65.  The  above  choice  of  the  kinetic  parameters  is  consistent 
with  the  cell  polarization  data  at  30-90  °C,  0.1-0.5  atm  P0l, 
35-100%  RH  (<2>)  and  0.054-0.186  mg  cm"2  Pt  loading  (Lpt)  and  also 
with  the  mass  activity  and  specific  activity  measurements  in  H2/O2  at 
0.9  V  IR-corrected  cell  voltage,  80°C,1  atmPH2  andPO2,andl00%RH. 
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Appendix  1.  Average  H2  and  O2  pressures  in  the  Cell 
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Fig.  12.  Sample  axial  profiles  of  y>a ,  PH]  and  <I>  in  run  1.2  at  0.4  A  cnr2. 


<,0  =  <,in(l  -  VSRc) 


<,0  =  <Jn(l  -  1/SRa) 


where  X  is  the  mole  fraction  and  the  subscripts  in  and  o  denote  the 
cell  inlet  and  outlet,  respectively. 


Appendix  2.  Representative  axial  profiles 

Fig.  12  presents  a  sample  of  the  calculated  axial  profiles  of 
i,Po2 ,  Ph2  and  #  for  Run  1.2  at  0.4  A.cm-2.  The  anode  H2  stream  in 
Fig.  12  flows  from  left  to  right  while  the  cathode  air  stream  flows 
from  right  to  left. 
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